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Abstract. Ground-state scalar isoscalar mesons and a scalar glueball are described in a U(3) x (7(3) chiral 
quark model of the Nambu-Jona-Lasinio (NJL) type with 't Hooft interaction. The latter interaction 
produces singlet-octet mixing in the scalar and pseudoscalar sectors. The glueball is introduced into the 
effective meson Lagrangian as a dilaton on the base of scale invariance. The mixing of the glueball with 
scalar isoscalar quarkonia and amplitudes of their decays into two pseudoscalar mesons are shown to be 
proportional to current quark masses, vanishing in the chiral limit. Mass spectra of the scalar mesons and 
the glueball and their main modes of strong decay are described. 

PACS. 12.39.Ki Relativistic quark model - 12.39.Mk Glueball and nonstandard multiquark/gluon states 
- 13.25.-k Hadronic decays of mesons - 14.40.-n Mesons 



1. Introduction 

Scalar mesons play an important role in the strong inter- 
actions of elementary particles . They are also signif- 
icant for a correct description of nuclear interactions || . 
Today, 19 scalar resonances are observed in the energy in- 
terval from 0.4 to 1.7 GeV However, their nature is not 
yet quite clear. Are they two-quark or four-quark states? 
Which of them belong to the ground state nonet of the 
U (3) flavour group; and which, to the excited one? Is there 
a glueball among them, and where does it lie? These are in- 
triguing questions, where to find answers, many attempts 
are currently being made by different authors. From a 
quick survey of what has been done by now, one could 
conclude that we do not yet approach the solution. A lot 
of modelsjsee, e.g., |,|0,|,|,00,|l],|l3|,|l4[|l|,|l|,0 
[l8|,|l9,20,|2l],^2|]) has been suggested, but none of the ap- 
proaches has given us the clue that could allow us to clar- 
ify the true nature of scalar mesons once for all. Moreover, 
one might be disappointed by the fact that different mod- 
els give different results that do not much overlap, being 
in some cases controversial. 

In our work, we describe scalar isoscalar mesons and 
their mixing with scalar glueballs. As we learn from QCD, 
if there were no quarks in the world, gluons themselves 
could form bound objects due to a large coupling strength 
of gluon self-interactions. This stimulated the search for 
bound gluon systems both in experiment and in theory. 
As the gluon carries no flavour, lepton or baryon quan- 
tum number, glueballs must be searched among isoscalar 
mesons. Indeed, the simplest gluonic formation possesses 
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the quantum numbers of a scalar isoscalar meson. If one 
looks into tables of experimental data, one finds a large 
number of scalar isoscalar states that can be mixed states 
of quarkonia, multi-quark systems, hybrids, and glueballs. 
To distinguish a glueball is really a difficult task, because 
we have no reliable test that would give us the truth. 

As the perturbative approach does not work here, dif- 
ferent phenomenological models and lattice simulations 
are involved in the study. From recent results |23|, f24| , |25f| 
one can conclude that it is most probably that glueballs 
are real objects of our world. There exist numerical esti- 
mates for probable masses of glueballs, however still in the 
world without quarks. Lattice calculations report that the 
lightest scalar glueball should be found between 1.5 and 
1.7 GeV. 

Amsler |ll| considered the state /o(1500) as a can- 
didate for the scalar glueball. QCD sum rules ]l9) and 
the K-matrix method [gfj showed that both /o(1500) and 
/o(1710) are mixed states with large admixture of the 
glueball component. Moreover, QCD sum rules |l9| re- 
quire that light glueballs (below 1 GeV) should exist, 
which is in contradiction with what lattice calculations 
suggest. 

A glueball cannot be searched without investigating 
the nature of the rest of scalar mesons that are not heav- 
ier than, 2 GeV and which we consider mostly as formed 
by quark-antiquark pairs. All the bound isoscalar qq states 
are allowed to mix with glueballs, and their spectrum has 
many interpretations made by different authors. For in- 
stance, Palano |^6| suggested a scenario, in which the 
states a (980), A'*(1430), / (980), and / (1400) form a 
nonet. The state /o(1500) is considered as the scalar glue- 
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ball. Tornqvist et al. [glj looked upon the states /o(980) 
and /o(1370) as manifestations of the ground and excited 
ss states; and the state /o(400 — 1200), as the ground 
uu state. Van Beveren et al. considered the states 
/o(400 — 1200) and /o(1370) as uu ground states; and 
the states /o(980) and /o(1500), as ss ground states. Two 
ground states for each qq system occur due to pole dou- 
bling, which takes place for scalar mesons in their model. 
Shakin et al. Jll| obtained from a nonlocal confinement 
model that the /o(980) resonance is the ground uu state, 
and /o(1370) is the ground ss state. The state /o(1500) is 
considered as a radial excitation of / (980). They believed 
the mass of the scalar glueball to be 1770 MeV. 

In our recent papers Jl7j , following the methods given 
in Refs. J|,|2^,^9],]3(|, we showed that all experimentally 
observed scalar meson states with masses in the inter- 
val from 0.4 to 1.71 GeV can be interpreted as members 
of two scalar meson nonets — the ground state of the 
meson-nonet and its first radial excitation. We considered 
all scalar mesons as qq states and took into account the 
sing let-octet mixing caused by the 't Hooft interaction. In 
|0, we obtained a scalar isoscalar state with mass 1600 
MeV and had to choose, to which of the experimentally 
observed states /o(1500) and /o(1710) we should ascribe 
it. From our analysis of the decay rates calculated in our 
model, we found that /o(1710) better fits to the nonet 
of quarkonia than /o(1500). Therefore, we supposed that 
the state /o(1500) contained a significant component of 
the scalar glueball (see ]l9[|o)). However, the final deci- 
sion should be made after including the scalar glueball into 
the model, and taking account of its mixing with quarko- 
nia, which will shift scalar meson masses. 

At present, there exist two candidates for the glueball: 
/o(1500) and / (1710) ||],|^||. To describe the proper- 
ties and mixing of the glueball with the other scalar states, 
one should introduce an additional scalar isoscalar dilaton 
field x m t° our model, in addition to the quarkonia which 
have already been described (l7). For this purpose, one 
can make use of the idea of approximate scale invariance 
of effective Lagrangians based on the dilaton model. Such 
models were studied by many authors (see, e. g., fl^,|22l,[3l"l, 
|32|,|33|]). Unfortunately, there is no unique way to introduce 
the dilaton field into a chiral Lagrangian. This justifies the 
large number of models dealing with glueballs. 

The guideline, one should follow when introducing the 
dilaton field into an effective meson Lagrangian, is to re- 
produce the Ward identity connected with the scale anomaly. 
The latter leads to the following equation for the vacuum 
expectation value of the divergence of the dilatation cur- 
rent 



Lagrangian by requiring that, in the chiral limit, our La- 
grangian should be scale-invariant except for the dilaton 
potential. To realize this program, one should multiply all 
dimensional parameters of the original Lagrangian (with- 
out dilaton) by a corresponding power of the dilaton field 
divided by its vacuum expectation value Xc to preserve 
the dimensions of model parameters. Thus, instead of the 
four-quark coupling constant G, the 't Hooft coupling con- 
stant K, ultraviolet cutoff A (necessary for the regulariza- 
tion of divergent integrals coming from quark loops), and 
the constituent quark masses m q (q = u,s), one should 
use G(xc/xf, K( X c/x) 5 , A(x/Xc) and m q (x/Xc)- 

Current quark masses are not multiplied by the 
dilaton field and violate scale invariance explicitly, as it 
takes place in QCD. Their contribution to the divergence 
of the dilatation current is determined by quark conden- 
sates and disappears in the chiral limit (see ([!])). 

Omitting, for a moment, the 't Hooft interaction in our 
approach (which, to an extent, is in the spirit of papers [ fTo[ 
^2|,^9| ), we require that the Lagrangian is scale-invariant 
in the chiral limit both before and after the spontaneous 
breaking of chiral symmetry (SBCS), except for the dila- 
ton potential. This property can be obtained by consid- 
ering (after bosonization when the effective Lagrangian 
is expressed in terms of bosonic scalar and pseudoscalar 
fields (7 and 6) the shift of the scalar meson field a 



a = a — rn- 



Xc 



(to = 0), 



(3) 



mJqq), 



(1) 



where (er')o = 0, (<r)n = —to, guaranteing that the rela- 
tion ([[]) is satisfied [||. The nonzero vacuum expectation 
value of a appears as a result of SBCS, and thus, the 
constituent quark mass is produced. In the case of non- 
vanishing current quark masses, (|^) changes by including 
an additional (non-scaled) mass term to into the r.h.s. 

This change produces an interaction term ~ (^-^ 

in the effective Lagrangian ( |3l| ) which breaks both chiral 
and scale symmetry just in the way required by the quark 
mass term m°qq of the QCD Lagrangian. 

The structure of the paper is as follows. In Section 2, 
we derive the usual U(3) x J7(3)-flavour symmetric effec- 
tive Lagrangian with the 't Hooft interaction and without 
dilaton fields. In Section 3, the dilaton field is introduced 
into the effective Lagrangian obtained in Section 2. Gap 
equations are investigated in Section 4. In Section 5, we 
derive mass terms and fix the model parameters. The main 
decays of scalar isoscalar mesons are calculated in Section 
6. Finally, in the Conclusion, we discuss the obtained re- 
sults. 
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Nr. 



1 



N 



-G 



24' " 12' J J W~'"7 1 ™ 
where N c is the number of colours; Nf, the number of 
flavours; ("Gj^) and (qq), the gluon and quark conden- 



sates; to , the current quark mass. 

In this paper, we are going to use the most natural 
method of introducing the dilaton field into the effective 



2. Chiral effective Lagrangian with 't Hooft 
interaction 

A U(3) x U(3) chiral Lagrangian with the 't Hooft inter- 
action was investigated in paper J34|. It consists of three 
terms as shown in formula (E). The first term represents 
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the free quark Lagrangian, the second is composed of four- 
quark vertices as in the NJL model, and the last one de- 
scribes the six-quark 't Hooft interaction J|5| that is nec- 
essary to solve the J7a(1) problem. 



G 

L = q{iB - m°)q + - $}(«A„g) a + (qij 5 \ a q) 2 } 

a=0 

-K{det[q(l + 75)9] + det[g(l - l5 )q]} . 



(4) 



Here G and K are coupling constants, A a (a = 1, ...,8) 
are the Gell-Mann matrices Ao = \/2/3 1, with 1 being 
the unit matrix; m is a current quark mass matrix with 
diagonal elements m°, m°, m° 8 (m° « m° d ). 

The standard bosonization procedure for local quark 
models consists in replacing the four-quark vertices by 
Yukawa couplings of quarks with bosonic fields which en- 
ables one to perform the integration over quark fields. The 
final effective bosonic Lagrangian appears then as a result 
of the calculation of the quark determinant. To realize 
this program, it is necessary, using the method described 
in p^, ^5[|3^ ,|37||, to go from Lagrangian ([|) to an interme- 
diate Lagrangian which contains only four-quark vertices 



L = q(id - m°)q + ~ £ [G^ (qr a q)(qnq) 

a,b=l 

+ G ifc } (075T o g) (qij5nq)], 



(5) 



where 



T a = A a (o=l, 7), T S = (V2X + Ag)/v3, 

r 9 = (-A + V2A 8 )/V3, 

G { n = Gg } = G ( 3 f =G±4Km s lf(m s ), 

G ( m = = G { 6 t ] = Gff =G± 4Km u lf(m u ), 

G$ = G t 4Kmj£{m s ), G^ = G, 

G ( i } = G$ = ±4V2Km u lf{m u ), 

G { a V=0 (^¥=b; o,6=l,. ..,7), (6) 



< = K - 32Km u m s I?(m u )I?(m s ), (7) 
m° = m° s - 32Km 2 J^(m u ) 2 . (8) 

Here m u and m s are constituent quark masses and the 
integrals 

InM=j^[dik ^f-^l (n = l,2; a=u,s), (9) 



where m a represents a corresponding constituent quark 
mass 1 : m u or m s . Note that we have introduced the no- 
tation of constituent quark mass already here, although 
they will be consistently considered only later, when dis- 
cussing mass gap equations (compare ( ffl| ) and (Q)) and 
the related shift of scalar meson fields (see ( p"2| ) and (p"3"|)). 
However, as we want to use an effective four-fermion in- 
teraction instead of the original six-quark one, we have to 
calculate quark loop corrections for the constant G (see 
(^|)) using full quark propagators with constituent quark 
masses. 

In addition to the one-loop corrections to the constant 
G at four-quark vertices, we have to modify the current 
quark masses mP a (see (fy and (Jsj) ) . This is to avoid the 
problem of double counting of the 't Hoot contribution 
in gap equations which was encountered by the author 
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After the redefinition of the constant G and of 
the current quark masses, we can guarantee that in the 
large- iV c limit the mass spectrum of mesons and the gap 
equations, derived from the new Lagrangian with modi- 
fied four-quark vertices, are the same as those which are 
obtained from the original Lagrangian with six-quark ver- 
tices. 

Now we can bosonize Lagrangian (^). By introducing 
auxiliary scalar a and pseudoscalar <f> fields, we obtain |3, 



CM) = -\ E (°a{G^)-J* h + UG {+) T a h 

a,b=l 

-1 Trln lid-m + ^T a {a a +i 1 * > 4> a ) \ (11) 



k 0=1 J 

As we expect, the chiral symmetry is spontaneously bro- 
ken due to the strong attraction of quarks in the scalar 
channel and the scalar isoscalar fields acquire nonzero vac- 
uum expectation values (a a )o 7^ (a = 8, 9). These values 
are related to basic model parameters G, m° and A via gap 
equations as it will be shown in the next Section. There- 
fore, we first have to shift the a fields by proper values so 
that the new fields have zero vacuum expectation values 



Ma 



Ma 1 



= 0, 



(12) 



where [i a = 0, (a = 1, . . . , 7), ^ 8 = m u , fxg = -ffl,/V2 
and M a = 0, (a = 1, . . . , 7), 0° = m° u , $ = -m°/V2. 
After this shift we obtain: 



(2tt)V e {k 2 +m 2 a ) n: 



C(a',<j))= L G {a',<j)) -i Trlnji<9- 
9 / 1 



(13) 



are calculated in the Euclidean metric and regularized by 
a simple 0(4)-symmetric ultraviolet cutoff A. For I^(m a ) 
one gets 



(10) 



where 



LgW, 



The notation "constituent" quark mass refers here to the 
total quark mass appearing in the full quark propagator. 
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a,b=l 



.,6=1 



From Lagrangian ( pL3[ ) we take only those terms (in mo- 
mentum space) which are linear, squared, cubic and quar- 
tic in scalar and pseudoscalar fields .0 

C(a',c/>) = L G (a',<j>) +tr[^(m)(a' 2 + 2 ) 

-AmI A {m)<x' + 2I A (m)(a' 2 + <j) 2 ) 
-Am 2 I A {m)cr' 2 + AmI A ( m y (a' 2 + </> 2 ) 2 
-I A (m)(a' 2 + 2 ) 2 + I A (m)[a' - m, 0] 2 _, (15) 



0-'=^CT Q T a , = ^ 0a T a , 



(16) 



a=l 



a=l 



where "tr" means calculating the trace over r-matrix ex- 
pressions and [•••]- stands for a commutator The ex- 
pression for I A (m a ) in Euclidean metric is given in (p"o[). 
The integrals I A (m a ) are also calculated in Euclidean space- 
time 



16tt 2 



Id I Ar + l I - 



T 2 



yl 2 



(17) 



Then, we renormalize the fields in (15) so that the kinetic 
terms of the effective Lagrangian are of conventional form, 
and diagonalize the isoscalar sector. 



£(a r ,<f> r ) = L G (a T 
+tr 



r„2 
P tr 2 



ir 2\ 



4mgI A (m)cr r 



+2g 2 I A (m)(a r 2 + Zcff 2 ) + -[to, ft 



-m 2 a r2 +mg<j r (a r2 



z<t> 



r 2\ 



[m,<l> r ]-[<T r 



-L((a r2 + Zclf 2 ) 2 -[a r ,clf} 2 _) 



a=l 



a=l 

For Lq we have: 

1 9 - 1 

~2 £ ida^a ~ + Ma) ( G(_) ) q 6 ^ 6 °* ~ M& + ^ 



,,6=1 
9 



2>«( G,+, L' 

1,6=1 



(20) 



2 Despite that the scalar fields are of the main interest in 
this paper, we still need pseudoscalar fields to fix the model 
parameters. 



Here we introduced Yukawa coupling constants g a - 

°'a = 9a<Ja, 4>a = VZg a (/> r a , (21) 



(18) 
(19) 



2 2 2 2 2 I a tA( M— 1 

01 = .92 = 5 3 =9&=9u = [4ia , 
04 = .9 2 = .9e = 5? = [^(rm,, m,)] -1 , 
ff9=52 = [^K)]- 1 , 



(22) 



I A (m u ,m s 



6»(yl 2 



(2tt) 4 J e (fc 2 + m 2 )(fc 2 +m 2 ) 



(47t) 2 (to 2 — to 2 ) 

2 i„ I i 2 

-l 



-m„ In | — 



1 



6to m 



to c In I — 7 

to; 



1.44, 



(23) 
(24) 



where we have taken into account Tr-^-transitions leading 
to an additional Z factor, with Ma 1 being the mass of 
axial- vector meson (see [21). The renormalized scalar and 
pseudoscalar fields in (121 ) are marked with the superscript 
r. 

The mass formulae for isovectors and isodublets fol- 
low immediately from (ff8|). One just has to look up for 
the coefficients at a r 2 and (ff 2 . There are still nondiag- 
onal terms in (^0|) in the isoscalar sector. This problem 
is solved by choosing the proper mixing angles both for 
the scalars and pseudoscalars (see e.g. As we are go- 
ing to introduce the glueball field, the mixing with scalar 
isoscalar quarkonia will change the situation. One has to 
consider the mixing among three states, which cannot be 
described by a single angle. For simplicity, in our estima- 
tions we resort to a numerical diagonalization procedure, 
not to the algebraic one. Concerning the pseudoscalar sec- 
tor, one can avail oneself with the results given in p4| . All 
what concerns dealing with the glueball is discussed in the 
next Section. 



3. Nambu-Jona-Lasinio model with dilaton. 

As we have already mentioned before, we introduce the 
glueball into our effective Lagrangian, obtained in the pre- 
vious Section, as a dilaton. For this purpose we use the 
following principle. Insofar as the QCD Lagrangian, in 
the chiral limit, is scale invariant, we suppose that our 
effective meson Lagrangian, motivated by QCD, has also 
to be scale invariant both before and after SBCS in the 
case when the current quark masses are equal to zero. 
Note that the scale anomaly of QCD is reproduced by 
the dilaton potential. As a result, we come to the follow- 
ing prescription: the dimensional model parameters G, A, 
and K are replaced by the following rule G — > G(xc/x) 2 i 
K — > -ftT(x c /x) 5 , A — ► A(x/Xc) 2 , where x is the dilaton 
field with the vacuum expectation value Xc- Moreover, 
the constituent quark masses are replaced by the rule 
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m a — ► m a (x/Xc)- Concerning the current quark masses, 
the are left unsealed. This leads to the following formula 

°a = o' a ~ (Ma ~ Ma + Ma) — + Ma 

Xc 



(Ma - Ma) — 

AC 



uP — 



(25) 



where /x° = 0, (a = 1, . . . , 7), (4 = m^s Mg = -m° s /V2, 
and the definition of M° is given after (U3). The difference 
Ma — Ma i s proportional to the 't Hooft interaction constant 
K (see (0) and (|J)) and has conventional scale behaviour, 
therefore it should be scaled in the same way as fi a - Note 
that in the r.h.s of ( ^5|) x' denotes the quantum fluctua- 
tions of the dilaton field around its vacuum expectation 
value Xc- 

Finally, we come to the following Lagrangian: 
C(a r , <j> r , x) = C{x) + L kin (a r , f) + L G (a r , <j> r , X ) 

3 



-tr 



-Amgl^{m)a r 



, X r / r 2 

+mg — er (a 

Xc 



Xc 
r 2\ 



X \ 2 2 r 2 

— —ma a 

Xc 



b r*)-2-{ a r2 

4 



2 f 



Here £(x) is the pure dilaton Lagrangian 
1 



£(x) 

with the potential 



-Ad vX y-v{x) 



V(x) = B ~ 
\Xo 



lnl^l -1 

Xo 



(26) 



(27) 



(28) 



that has a minimum at x = XOi an d the parameter B 
represents the vacuum energy, when there are no quarks. 
The curvature of the potential at its minimum determines 
the bare glueball mass 



4-/B 
Xo 



(29) 



The part L\ tin {a r , <jf) of Lagrangian (|26|) contains pure follows 
kinetic terms 



(30) 



2 9 



VAC/ o,6=l 



1 fx 



2 \X 



2 9 



,6=1 



Ox X 



u° — 



(g(- 



-1 

ab 



/ r I -0\ X o X 



2 9 



(31) 



a, 6=1 



The dilaton field is here expanded around its vacuum ex- 
pectation value: x = x' + Xc (x)o = Xc (x')o = 0. 

Recall that the terms proportional to m break ex- 
plicitly chiral and scale invariance in the same way as 
the current mass term of QCD Lagrangian. Notice also 
that for our linear cr-model (|2^), together with the gap 
equation (see (ffl|)), lead to a scale-invariant pion term 
~ M%ir 2 (x/Xc) /2 instead of the scale-violating term 
~ M^7T 2 (x/x c )/2 arising in nonlinear cr-models |}3,38 . 

As one can see, expanding (|3l]) in power series of X , 
we can extract a term that is of order x 4 - It can be ab- 
sorbed by the term in the pure dilaton potential which 
has the same degree of X . Obviously, this leads only to 
a redefinition of the constants B and X o-> which anyway 
are not known from the very beginning. Moreover, saying 
in advance, terms like % 4 do not contribute to the diver- 
gence of the dilatation current ([!]) because of their scale 
invariance. 

Let us now consider the vacuum expectation value of 
the divergence of the dilatation current calculated from 
the potential of the effective meson-dilaton Lagrangian 



AB 




E m ?(w)- 

q— u,d,s 



a' —0 



(32) 



Here V = V( X ) + V (a r , <f) r , x) , where V \a r , (jf , X ) is the 
potential part of Lagrangian £(a r , (ff , X ). Note that we 
have simplified (^), taking into account that the quark 
condensates are related to integrals Ii(m u ) and I\{m s ) as 



4m q Ii(m q ) = -(qq) , (q = u,d,s), 



(33) 



and that these integrals are connected with constants G 



(-) 



to which we pay no further attention. The next term reads through gap equations, as it will be shown in the next 

Section. Comparing the QCD expression ([!]) with (|32"|), 

L G (a r ,<f> r , X ) = " ' ~ 

Ef^<-(^-^ + Ma)-+M° , l (G { - 
i A Xc 

a, 0=1 



one can see that the quark condensates enter into both 
formulae in the same way. Equating the right hand sides 
"of (HI) and (IF 



a b 



x ( 9wl ~ (Mb - Mb + Mb)— + Mb 

V Xc 



Cg~ E rn° q (qq)=4B 

q—u,d,s 



E m °^> ( 34 ) 

q—u,d,s 
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we obtain the correspondence 



C g = AB 



(35) 



This equation relates the gluon condensate, whose value 
we take from other models (see e.g. j39|), to the model 
parameter B. The next step is to investigate the gap equa- 
tions. 



4. Gap equations 

As usual, gap equations are follow from the requirement 
that the terms linear in a r and x' should be absent in our 
Lagrangian 



SC 



,x')=o 



SC. 



SC 
= 0. 



(36) 



<r r ,x')=o 

This leads to the following equations 



m°)(G("))ri 



-(G ( - } ) 



Nil 



V2 

-8m u I?(m u )=0 ,(37) 
(m s - m° s )(G^ - V2(m u - fhP u )(G^i 

-8m s Ii(m s ) = 

4 



AB ^ 



-iln 

Xo 



Xc 



,(38) 
(39) 



Here 



A = - 
2 



9 

£ 

a,b=l 



(40) 



and 



is proportional to the current quark masses fi® ~ m b 
thereby small. 

Using (j^) and (^J), one can rewrite the gap equations 
@ and (||) in a well-known form |37f 



m° = m u - 8Gm„J^(rn„) 

-32ifTO„m s // 1 (TO u )/^(m s 
= m s — 8Gm s I^(m s ) 
-32K(m u lf{m u )) 2 . 



(41) 
(42) 



The equations discussed above allow us to relate the 
current quark masses to the rest of model parameters and 
also to relate the constants B and xo to the gluon con- 
densate and Xc- The constituent quark masses, ultraviolet 
cutoff, and four-quark coupling constants will be fixed, as 
usual in NJL, by means of the Goldberger-Trcimann re- 
lation, the p — ► 7T7T decay constant, pion weak decay con- 
stant and the mass spectrum of pseudoscalars (For details 
see |34| and Refs. therein). In the next Section we define 
Xc, using the bare glueball mass (without mixing effects) 
as a parameter. 



Table 1. The masses of physical the scalar meson states 07, 
on, on 1 and the values of the parameters Xc, Xo, bag constant 
B, and (bare) glueball mass m g (in MeV) for two cases: 1) 



1500 MeV and 2) M aiII = 1710 MeV. 





01 


on 


OH l 




Xo 


B, [GeV 4 ] 


m 9 


I 


555 


1075 


1500 


191 


192 


0.005 


1480 


II 


555 


1080 


1710 


167 


168 


0.005 


1695 



5. Mass formulae and numerical estimations. 

The potential part of Lagrangian ( p6| ) which is quadratic 
in fields a r and x' an d which we denote as has the 
form 



l( 2 V,^,x') 

1 



T.9n 



{[(G<- 
{[(G<- 



J 88 



sift™,,)] 



r2 



4m 2 > 



J 89 cr 8 cr 9 



2 C 



a,b=8,9 A 



)ab9bVbX' 



(43) 



The dilaton and its interaction with quarkonia does not 
change the model parameters m u , m s , A, G, and K fixed 
in our earlier paper |34[ 

to„ = 280 MeV, m s = 420 MeV, A = 1.25 GeV, 
G = 4.38 GeV" 2 , K = 11.2 GeV -5 . (44) 

As it has been already mentioned, after the dilaton field is 
introduced into our model, there appear three new param- 
eters: xo, Xc, and B. To determine these parameters, we 
use the two equations ( |35| ) and ( |39| ) and the bare (without 
mixing effects) glueball mass 



AC a -A 



(45) 



We adjust it so that, in the output, the mass of the heavi- 
est meson would be 1500 MeV or 1710 MeV, and thereby 
fix Xc - For the gluon condensate, we use the value (390 MeV) 4 
[ p9[ . The result of our fit is presented in Table [l] where we 
show the spectrum of three physical scalar isoscalar states 
<ji, a 1 1 and am. The last one is associated with the glue- 
ball. The parameters xo and B are fixed by the gluon 
condensate and constituent quark masses 



Xo = Xc exp 



B 



• exp 



A 

2C~ g 
2A 



(46) 
(47) 



The mixing of scalar isoscalar fields is described by the 
matrix b that connects the nondiagonalizcd fields o r — 
(<7g,<Tg,x') with the physical ones a p h = (pi,on,oni) 

o r =ba ph . (48) 

The matrix elements of b are given in Table ||. 
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Table 2. Elements of the matrix b, describing mixing in the 
scalar isoscalar sector. The upper table refers to the case (Tin = 
/o(1500), the lower one to the case am = /o(1710) 





<TI 


a 1 1 


(Till 




0.9804 


0.1865 


-0.0636 




-0.1963 


0.9535 


-0.2288 


1 

X 


0.0180 


0.2368 


0.9714 






an 


(Till 


&u 


0.9804 


0.1912 


-0.0474 




-0.1965 


0.9672 


-0.1609 


f 

X 


0.0151 


0.1671 


0.9858 



large, compared to the pion case as wP s m°. In the same 
way, the amplitude can be split into two contributions 



A a UI ^KK - A l Ilt ->KK 



A q 



(54) 



where the pure glueball decay into KK is represented by 
amplitude 



A 9 



2M| 



(55) 



Its value is large and comparable with the quark compo- 
nent contribution 



6. Decay widths 

Once all parameters are fixed, we can estimate the decay 
widths for the main strong decay modes of scalar mesons: 
u\ — ► 7T7T, cr / — ► KK, a i — > rjrj, ai — * 7777', and er; — -> Air 
where I = I, II, III. 

The amplitudes that describe the decays are relatively 
simple in our model. The decays of quarkonia were con- 
sidered in Here we only give numerical estimates 
for their decay widths, where the mixing with glueball is 
taken into account (see Table |2|). Below, we discuss only 
those amplitudes that describe glueball decays. The pro- 
cess (J in — > tttt is given by the amplitude 



A —A3 

■"■"ill — *mt ■ rL ain—>irir 



A q 
am — ^7T7r 



(49) 



which has been divided into two parts. The first part rep- 
resents the contribution from the pure glueball. It is pro- 
portional to the square of the pion mass 



M 2 

A3 - —E h 

"■a 1 1 j— urn u X<?iii- 
Ac 



(50) 



where b x<7lII represents a corresponding element of the 
3x3 mixing matrix for scalar isoscalar states (see Table ||). 
This contribution is small (since it is proportional to the 
current quark mass m°), and the process is determined 
by the second part that describes the decay of the quark 
component of the glueball 



A q 

(T/j/— urir 



= 2g u m u Zb 



(51) 

<C 1, this term 



Despite the smallness of mixing, \b au , 
prevails over the pure glueball contribution because M 2 
is noticeably less than 2g u m u Zx c ba u cr III ■ As a result, the 
decay width of om, if is 



A q 

trin—>KK 



2g u m u Zb (Tu 



2V2<7, 



s m s Zb asCJ[II . 



(56) 



In this case, the contribution from the quark component 
is provided by both u{d) and s quarks. In the case that 
&IIJ is /o(1500), we have 



r aiII ^KK = & MeV, 
and in the other case (07 // = /o(1710)) 
r aui ^ Kit = 90 MeV. 



(57) 



(58) 



Strange quarks contribute more and interfere with the 
pure glueball part, essentially reducing the decay width 
(by a factor 3). 

The amplitude of the decay of a glueball into 7777 and 
7777' can also be considered in the same manner. The only 
complication is the singlet-octet mixing in the pseudoscalar 
sector. The corresponding amplitudes are 



A ctti — »m — A o-in^vr) 



A q 

triii—njri 



A 9 



M 2 

Xc * 



(59) 
(60) 



A liu^vv = 2 9u m uZb,y u!yiII sin 8 

-2V2g s m s Zb asaiII cos 2 8, (61) 

where 8 = 6 — 80, with 8 being the singlet-octet mixing 
angle in the pseudoscalar channel, 8 w —19° Q, and 8q 
the ideal mixing angle tan#o = \j\pl. The decay widths 
thereby are: 

r„ lu ^ m = 25 MeV 
for am = /o(1500), and 



r„ 



42 MeV 



(62) 
(63) 



r ain ^ = 4 MeV 
for g in = / (1500), and 



(52) for Gin = /o(1710). 

For the decay of the glueball into 7777', we have the 
following amplitude 



, ff7r = 3 MeV 



(53) 



-4, 



triii— >rfr)' 



A 9 an ^ w + Al u ^ vv „ (64) 



for Gin = /o(1710). As one can see this process occurs 
with a relatively low rate. 

In the case of KK channels, the contribution of the 
pure glueball is also proportional to the mass square of 
the secondary particle, kaon in this case. But it is rather 



A 9 

<Tin—*riri' 



0, (65) 

—2Z sin 2d(g u m u b au a IIJ 

+V2g s m s b (TBCriII ). (66) 
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The amplitude A 9 aiii , is equal to zero because there is 
no decay of a bare glueball into 7777'. This process occurs 
only due to the mixing between the glueball and scalar 
isoscalar quarkonia. The decay widths are as follows 



r, 



fill— ►TJJj' 



for am = /o(1500), 



r, 



= 5 MeV, 



= 5 MeV. 



(67) 



(68) 



for am = /o(1710). The estimate for the decay /q(1500) 
into 7777' is just qualitative because the decay is allowed 
only due to the finite width of the resonance as its mass lies 
a little bit below the 7777' threshold. The calculation is made 
for the mass of / (1500) plus its half-width. For / (1710), 
we have a more reliable estimation since the mass is large 
enough for the decay to be possible. One can see that the 
order of magnitude for this decay is about 5 MeV. The 
estimate for /o(1500) is not in contradiction with it. 

The decays into four pions are estimated as decays pro- 
ceeding throug two channels: one with two intermediate 
scalar resonances (x — * era — > Air) and one with only one 
intermediate scalar resonance (x — > a2ir — > <tir). Here we 
neglect the mixing of glueball with quarkonia since the 
mixing effect is small (it is proportional to 777,").^] The am- 
plitude describing the decay into 2tt + 2it~ is as follows 



.4 



(7 1 1 1 — *27T' 2 "71 



^■trj//^(T(T^27r+27r- + Arm -^(t2tt^2tt+2ti- j 



(69) 



.4 



2>2mtZMl 



<7 j j j — *tr (J — >27T~^ 2tt 



(A(s l2 )A(s u ) 



A 



(Tin — XT2-K— +2-7T+2-/T 



+A(s u )A(s 23 )), (70) 
16mf.Z . . . . . . . 

— £^ f -(^12) + ^(«34) 

+A(s u )+A(s 23 )), (71) 



where = 93 MeV is the pion week decay constant, Mo- 
is the mass of the state 07. The function A(s) appears due 
to the resonant structure of the processes 



(72) 



where F a is the decay width of the 07 resonance (see be- 
low) . This function depends on an invariant mass squared 



Sij defined as follows 



2 

3 ) 1 



(i,i = l,...,4). (73) 



Here i and j enumerate the momenta k, of pions n + (ki), 
ir~ (k 2 ), ir + (k 3 ), and 7r~(fc4). The amplitude describing the 
decay into 2it it + tc~ has the form 

— >27r (, 7r+7r~ 



3 Only the lowest scalar isoscalar resonance is taken into 
account here. The contribution from /o(980) should be notice- 
ably smaller because of a large mass and a narrow width of 
/o(980). 



-4 



(TJII — >(T(T — ►27r°7r+7T - 



.4 



A(s 12 )A(s 3i ),{75) 

16mf.Z 



(Jin — >{j2tt — >27r°7r+7r 



Xc 



(A(s 12 ) + A(s 34 )%76) 



Here k\ and k 2 are momenta of the two 7r°, and s\ 2 is their 
invariant mass squared. The indices 3 and 4 stand for 7r + 
and 7T _ , respectively. 

In the case of the decay into 47r° , we have 



A 



a 1 1 1 



A 



CF 1 1 1 — >(T(T — *-47T U 



.4 



Aa UI ^aa^4TV° + A aj u -^ a2 T:^A-n a > (77) 

16m 4 ZM 2 

+A(s l3 )A(s 24 ) + A(s 1A )A(s 23 )), (78) 
16mtZ 



a in -^xt2tt — »-47r 



+A(s 23 ) + A(s 24 ) 



A(s 13 ) + A(s u ) 
^A(s 34 )). (79) 



From our estimation it follows that in the case, where 
(T777 is identified with /o(1500), we have the total width 



Frj 



ilT = 30 MeV, 



and in the other case (am = /o(1710)) 
F aill ^ - 60 MeV. 



(80) 



(81) 



Let us present the decay widths of 07 and 07/. The 
state a 1 that we identify with /o (400— 1200) decays mostly 
into a pair of pions, and this process determines the width 
of 07: 

r ai ^„ « 760 MeV. (82) 

The state 07// does not affect it noticeably, since the mix- 
ing of the glueball with uu is very small. Therefore, the 
decay rate for both CT777 = / (1500) and 07// = /o(1710) 
is approximately the same in magnitude. 

The decay of the state an that we identify with /o(980) 
into pions is determined by the quark component and 
is slightly reduced by the glueball component because of 
mixing with the ss quarkonium. We obtain 



r an ^ = 17 MeV, 



if am = /o(1500) and 



15 MeV, 



(83) 



(84) 



if am = /o(1710). From experiment, we know that its 
decay width lies within the interval from 40 MeV to 100 
MeV. Concerning the process an — > nn, we obtain a decay 
width that is lower than the experimental one. Notice that 
this prediction is completely based on singlet-octet mixing 
following from the 't Hooft interaction j54j where dilaton 
effects do play a minor role. The decay into KK can also 
be taken into account. From experiment we learn that the 
decay into KK can contribute about 30% to the total 
width jj] . Our estimates for decays of the glueball are 
collected in Table |[ 
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Table 3. The partial and total decay widths (in MeV) of the 
glueball for two cases: am = /o(1500) and am = /o(1710), 
and experimental values of decay widths of /o(1500) and 
/o(1710) §. 





*7T7T 




r vv 


- 1 7777' 




-Tiot 


1 tot 


/o(1500) 


4 


42 


25 


5 


30 


100 


112 


/o(1710) 


3 


90 


42 


5 


60 


200 


130 



7. Conclusion 

As it was mentioned in the Introduction, the inclusion of a 
scalar glueball into the effective meson Lagrangian is quite 
an ambiguous procedure. The goal of our paper is to find 
the most physically justified way to do this. In the ap- 
proach presented above, we assume that (with the excep- 
tion of the dilaton potential) scale invariance holds for the 
effective Lagrangian before and after SBCS in the chiral 
limit. The terms depending on current quark masses break 
both the chiral and scale invariance, in accordance with 
QCD. This leads to the requirement that we should intro- 
duce the dilaton field into the constituent quark masses 
while the current quark masses remain unsealed. 

In this version of a scaled NJL model, the terms that 
describe mixing of the glueball with quarkonia are also 
proportional to current quark masses. The same is true for 
the amplitudes describing decays of the glueball into pairs 
of pseudoscalars. Insofar as the masses of current quarks 
are small in comparison with the other model parameters 
(constituent quark mass, \ c , A, and so on), this results in 
a small mixing of the glueball with quarkonia, relatively 
small rates for decays of the glueball into rnr, and only 
slightly changes the decay width of /o(980) — > 7T7T calcu- 
lated before introducing the glueball p3 . The decay of the 
glueball into two pions is mostly determined by its qq ad- 
mixture despite the small mixing. The mixing coefficient 
here, although being small (~ —0.06 if am = /o(1500)), 
is multiplied by a relatively large constant describing the 
decay of the cr-meson into a pair of pions. 

In the case of the KK channel, both the gluonic and 
quark components play an important role since the inter- 
ference between the gluonic and quark amplitudes is large. 
The relatively small contribution from the uu component 
slightly increases the decay rate of the glueball. But the 
contribution from the ss component reduces the contribu- 
tion from the pure glueball by factor 3. 

The decay into r/rj is mostly determined by the glueball 
component. The mixing of the glueball with ss reduces 
the decay rate but not significantly. The decay into 7777' is 
less than into 7777 and is allowed only due to the mixing 
of quarkonia with the glueball. This process serves as a 
measure of this mixing. However, in the case of /o(1500), 
it is difficult to give reliable estimates for its rate because 
the process occurs near the threshold. 

Decays into 4 pions are represented by two processes. 
In the first one, two intermediate scalar resonances are 
born by the glueball with their subsequent decay into two 
pairs of pions. In the second process, only one intermediate 
scalar resonance together with a pair of pions are produced 



immediately after the decay of the glueball. Then, the 
scalar resonance decays into pions. From our calculations 
it follows that the second process is dominant and two 
scalar resonances are less probable to appear. 

The total width of the third scalar isoscalar state is 
estimated to be about 100 MeV for M aiu = 1500 MeV 
and 200 MeV for M aill = 1710 MeV. If we assume that 
the /o(1500) state is the scalar glueball, the total decay 
width derived from our model is close to the experimental 
value (112 MeV). Unfortunately, the detailed data on the 
branching ratios of /q(1500) are not reliable and contro- 
versial |y. 

In conclusion, we would like to note that, in our model, 
the width of the decay of a glueball into two pions is small, 
because the amplitude describing this decay, is propor- 
tional to the current mass of it-quark (~ M% ~ m°) and 
does not depend on momenta. The latter in the chiral 
limit formally disagrees with the low-energy theorems ob- 
tained in paper |40(| . In general, we could consider a version 
of our model containing momentum-dependent vertices, 
whose momentum dependence is in agreement with these 
low-energy theorems. However, such a momentum depen- 
dence of the amplitude leads to too large decay width of a 
heavy glueball (see fl33|l), which contradicts the experimen- 
tal data. This witnesses to the fact that these low-energy 
theorems are not justified to be applied in the case of a 
heavy glueball. 

The results obtained here correspond to the leading 
order in 1/N C expansion (Hartree-Fock approximation). 
Next-to-leading order corrections can to an extent change 
the final results. Note also that, in the energy region un- 
der consideration (~ 1500 MeV), we work on the brim 
of the validity of exploiting the chiral symmetry that was 
used to construct our effective Lagrangian. Thus, we can 
consider our results as rather qualitative. Nevertheless, a 
satisfactory agreement with experimental data is obtained 
for the total width of / (1500). 

We are going to use this approach in our future work 
for describing both glueballs and ground and radially ex- 
cited scalar meson nonets which lie it the energy inter- 
val from 0.4 to 1.71 GeV. Small mixing angles make us 
hope that introducing the glueball into our model will not 
change the whole picture dramatically. 

We are grateful to Prof. S.B. Gerasimov, Dr. A.E. Dorokhov, 
and Dr. N.I. Kochelev for useful discussions. The work is sup- 
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program 2000 and the Graduiertenkolleg "Elementarteilchen- 
physik" of the Humboldt University, Berlin. 
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